eT BA Mee ht BY S Be Be) A Eat 


究 


万 美 言 ， 鲁 同 所 “， 雇 候 含 “， 于 和 白雪， 高 贝 贝 ， 


1. ”西藏 大 学 理学 院 ， 拉萨 850000 ; 2. 中 国 科 学 院 上 海 应 用 物理 研究 所 ， 上海 嘉定 201800 ; 


摘 要 : 本 文 利用 最 大 业 谱 估计 的 方法 研究 了 风云 二 号 D 星 >2.0MeyV 的 高 能 电子 通 量 的 周期 
特性 ， 该 算法 以 自 相关 (AR ) REKREI  WIMRAMEDS EHO , 推测 出 风云 二 

星 高 能 电子 通 量 存 在 13.87d 和 27.8d 的 周期 。 根 据 FPE、AIC 准则 确定 最 佳 阶 数 进而 计算 AR 模 

型 参数 ， 并 将 Levinson-Durbin、Burg BASRA MA AH THR |, AIRA ETE 
在 周期 特性 研究 上 存在 优势 。 这 一 结果 对 研究 地 球 同步 轨道 高 能 电子 空间 分 布 、 预 报 高 能 电子 增强 事 
件 和 预警 深层 充电 事件 至 关 重 要 。 


关键 词 : 高 能 电子 ; RAM: 周期 特性 ; 功率 谱 估计 
中 图 分 类 号 : P172.4 ”文献 标识 码 : A 文章 编号 xxxx - XXKX(KKKX)XK - X00 - XX 


高 能 带 ~ 组 成 的 空间 辐射 环境 是 影响 航天 器 在 轨 安 全 飞行 最 重要 的 空间 环境 之 一 ， 
主要 包含 ， 其 一 ， 这 些 高 能 带电 粒子 入 射 到 物质 中 ， 在 航天 材料 表面 或 内 部 发 生 充 放电 反应 ， 导 致 航天 
es 其 二 ， !， 改 变 元 器 件 内 部 的 电荷 分 
布 ， 造 成 元 器 件 的 逻辑 状态 混乱 而 引发 故障 甚至 失效 ， 其 三 ， 高 能 带电 粒子 的 辐射 效应 A 
作用 导致 航天 器 的 材料 、 元 器 件 性 能 变 差 ， 直 至 失效 ， 同 时 辐射 效应 也 会 对 人 体 造 成 损伤 ， 严 重 影响 
航天 员 的 安全 0。 

风云 二 号 是 我 国 研制 的 第 一 代 地 球 同步 轨道 气象 卫星 ， 与 极地 轨道 气象 卫星 相辅相成 ， 
构成 我 国 气象 卫星 应 用 体系 中。 风云 二 号 D 星 于 2006 年 12 月 8 日 成 功 发 射 升 空 ， 在 高 度 为 35783- 
35788km 的 高 轨道 上 飞行 ， 包 含 高 能 质子 、 电 子 、He 离子 等 ， 共 有 7 个 能 道 。FY-2D 定位 东 
经 86.5” 赤 道上 空 ， 轨 道 周期 为 1436.04 分 钟 ， 轨 道 倾 角 为 1.84 ， 自 旋 稳定 ， 转 速 100r/min。 风 
云 二 号 D 星 搭载 空间 粒子 探测 器 ， 它 是 空间 环境 监测 器 的 一 部 分 。 风 云 二 号 D 星 与 之 前 发 射 的 A 
k 同 实现 双星 观测 ， 发 挥 1+1 大 于 2 的 系统 效益 ， 提 高 卫星 观测 频次 ， 扩大 卫星 观测 范围 ， 
实现 立体 观测 和 动态 监测 叫 。 对 研究 高 能 粒子 的 周 期 特性 提供 更 精确 、 可 靠 的 数据 ， 同 时 也 可 
为 研究 太阳 活动 、 地 磁 活 动 对 空间 环境 的 影响 提供 很 好 的 参考 。 

本 文 主要 利用 最 大 炉 谱 估计 研究 了 近 地 空 间 高 能 电子 ( 宇 2.0MeV 电子 通 量 ) 的 周期 性 变 
化 特征 ， 同 时 对 比 其 与 Levinson-Durbin 和 Burg 算 法 的 优 缺 性 。 本 课题 研究 所 采用 的 数据 覆盖 
时 间 范 围 为 2006 年 12 月 19 日 至 2012 年 5 月 21 日 ， 时 间 精 度 为 Smin， 对 数据 空缺 部 分 均 进 行 插 
值 处 理 。 


RI 


* 收 稿 日 期 : xxxx 年 xx 月 xx 日 

资助 项 目 : 国家 自然 科学 基金 (11803024，11747128) , 西藏 自治 区 自然 科学 基金 项 目 (XZ2019ZRG-163) , 西藏 大 学 教师 学 历 提升 
“博士 后 计划 ”; 西藏 大 学 理学 院 2019 年 度 和 2020 度 院 级 教改 项 目 资助 

* 作者 简介 : 鲁 同 所 (1987-) , 男 ， 河 南 濮阳 ,西藏 大 学 理学 院 物理 系 ， 硕 士 生 导师 ， 中 国 科学 院 天 体 物 理学 博士 Email: 


1 


202011.00015v1 


chinaXiv 


lutongsuo@163.com 


ChinaXiv 合 作 期 刊 


chinaXiv:202011.00015v1 


ChinaXiv 合 作 期 刊 


1.1 mA AT 

NEDI TE Ee, TRO, AR CE FS 
性 上 傅 里 叶 变 化 缺乏 时 间 和 频率 的 定位 功能 ， 傅 里 叶 变 化 是 一 种 整体 变化 ， 它 不 能 够 反应 某 个 
时 间 段 的 信号 发 生 了 怎样 的 变化 ， 也 无 法 获得 某 一 频率 出 现 的 时 刻 信息 ; 其 二 傅 里 叶 变 化 对 非 
平稳 信号 存在 局 限 性 ， 本 文 研究 的 高 能 电子 通 量 ， 它 的 频率 是 随时 间 变 化 的 ， 用 傅 里 叶 变 换 进 
行 分 析 只 能 给 出 高 能 电子 通 量 频率 变化 的 整体 效果 而 不 能 完整 地 反映 出 某 一 时 刻 信号 的 本 质 特 
性 。 

最 大 入 谱 估计 具有 高 分 辨 率 和 短 时 性 的 特点 ， 其 原理 可 概述 为 :， 利用 已 有 的 自 相关 函数 
值 , WACK A BEE, 利用 N 个 已 知 的 自 相 关 函 数值 来 外 推 其 他 未 知 的 自 相 关 函 数值 ， 最 后 作 频 
域 变 换 , 求 得 连续 的 功率 谱 估 计 ， 为 研究 数据 的 周期 特性 提供 了 算法 文 撑 。 
AIRY”: 


计算 初始 值 
1 N=1 
"OO)= 二 > [xn] a) 
N To 


f(n) = 8 (n) = xn) 


求 出 预测 均 方 误差 pp 递 推 公式 


pz (1 bel) Pi (2) 


求解 AR 模型 的 反射 系数 kn 


N-1 
py faigna 一 了 
ke nm ,m=1 (3) 


Dfa DF + nD | 
求 前 向 预测 误差 和 后 向 预测 误差 ， 然 后 顺 次 估计 出 反射 系数 Kn 


n=m 


FEM) =f, (0) +k Sp. (nl) (4) 


Em (n) =f, a(n =1) T Eman -1) 
FH Levinsion 递 推 ， 求 出 阶 次 Mm = 2 时 AR 模型 参数 


AyD) = An {D+ Kim m1) i=1,2, =+, ml (5) 


m 


d,, M) = kn 
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重复 上 述 过 程 ， 直 到 m 等 于 所 需 AR 模型 阶 数 ， 求 出 所 有 的 AR 模型 参数 ， 然 后 预测 出 
WFAA R, B, 


m 


时 


; (6) 


RP k=1, 2,… ,Pp 为 p 阶 线性 预测 滤波 器 系数 ，c 为 预测 滤波 器 的 预测 误差 功率 。 


1.2 Levinson-Durbin 和 Burg 算法 
Levinson-Durbin、Burgts 算 法 均 是 线性 预测 方法 ， 是 根据 已 知 的 时 间 信 号 序列 计算 功率 谱 
估计 值 的 递 推算 法 。Burg 算法 利用 前 向 滤波 误差 有 ,和 后 向 滤波 误差 b, ,来 求 出 滤波 误差 功率 
为 最 小 的 a,，， 青 根据 Levinson-Durbin 算法 计算 AR 模 型 参数 aq, 。 要 利用 AR 模 型 进行 功率 谱 估 
计 ， 需 要 用 到 Yule-Walker 方程 : 
Gn i) = an Ò + kpaa (m-i), i= 1, 2,....m—l1 


M 阶 AR 模 型 的 第 m+1 个 参数 G 满 足 : 


G=p (7) 


pm 是 预测 功率 误差 ， 可 得 递 推 公式 : 


Py, = Paa- k, (8) 


由 〈8) 式 进行 递 推 ， 可 求 得 最 终 表 征 随机 信号 的 AR 模型 的 p+ 参数 ， 最 后 求 得 随机 信号 
的 功率 谱 密度 : 


Se”) = 0° o|H(e”)) 


1.3 AR( 自 相关 ) 模型 阶 数 的 确定 

AR 模 阶 数 的 确定 十 分 重要 ， 在 最 大 烂 谱 估 计 和 Levinson-Durbin, Burg 算法 中 都 会 用 到 。 
若 选 择 太 小 的 阶 数 ， 得 到 的 功率 谱 图 会 太平 滑 ， 不 能 够 有 效 地 分 辨 出 时 间 序 列 周 期 分 量 ， 阁 选 
择 太 大 的 阶 数 ， 会 影响 最 大 科 估计 值 的 稳定 性 。AR 模 型 的 最 佳 阶 数值 需要 在 递 推 过 程 中 来 确 
定 。 使 用 Levinson 弟 推算 法 ， 由 低 阶 到 高 阶 的 每 一 组 参数 及 模型 的 最 小 预测 误差 功率 pi, 是 
递减 的 。 理 论 上 来 说 ， 当 z 值 达到 期 望 值 ， 或 不 再 发 生 任何 变化 时 ， 即 满足 最 佳 适用 性 ， 此 
时 的 阶 数 就 是 应 选 的 最 佳 参 数 。 以 下 是 两 个 常用 的 阶 数 适 用 性 检验 准则 口 ; 
(一 ) 最 终 预 测 误差 准则 (FPE) 

FPE 准则 考虑 到 的 是 预测 值 和 真实 值 之 间 的 误差 ， 以 误差 最 小 为 原则 ， 而 使 误差 最 小 只 对 
应 一 个 阶 数值 > 由 此 AE r) 定义 为 : 


刊 
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map (9) 
N-(r+) 
式 中 N 为 高 能 电子 通 量 的 数据 样本 数 。 
(=) 赤 池 信息 准则 (ATO) 
AIC 准则 是 建立 在 烂 的 概念 上 ， 可 以 判断 所 估计 模型 的 复杂 度 和 此 模型 拟 合 数据 的 优良 
性 ， 定 义 AIC(r) Ai: 
AIC(r)=2r-21n(L)r (10) 


Bir SCA EAC AE RAG PPE (27) 和 AIC(2) 的 值 也 将 发 生变 化 。 在 z 取 菜 值 时 , PECON ATC (r) 
的 值 都 为 最 小 值 ， 那 么 此 时 所 取 的 7 就 是 最 佳 阶 数 值 。 


在 实际 运算 中 发 现 ， 数 据 较 短 时 ， 
阶 数 值 基本 一 致 。 即 ; 


所 得 的 最 佳 阶 次 值 偏 低 ， 且 以 上 两 个 准则 所 得 的 最 佳 


lim!g FPE() = AIC(7) (1 


No% 


2. 高 能 电子 通 量 周期 特性 


2. 1 数据 选取 


本 文选 择 分 析 能 道 宇 2. 0MeV 的 高 能 电子 通 量 ， 是 能 量 极 高 的 带电 粒子 。 由 于 初始 数据 在 
12 月 30 日 、2011 年 2 月 2 日 -2 月 11 日 、2011 年 8 月 22 日 -8 月 24 日 等 时 间 段 均 存在 
BRAG, 为 了 避免 缺失 数据 对 研究 周期 特性 的 影响 ， 因 此 利用 函数 Fillmissing 通 过 matlab 实 现 对 
所 收集 到 的 数据 进行 插值 处 理 ， 如 图 1 所 示 , 是 2006 年 12 月 19 日 到 2012 年 5 月 21 日 高 能 电子 在 能 


2009 年 12 月 6 日 


1) 


WH S2. 0MeV 的 电子 通 量变 化 情况 , 数据 频次 为 5min, 每 日 24 小 时 ， 平 均 每 日 288 个 数据 ， 经 过 插 


值 处 理 数据 后 ， 数 据 量 超过 50 万 。 


4 
a, 


wo 
T 


22.0 MeV Electron Flux 
(cm? s-sr-keV) 1 
= N 


0 


DO e] 
N ss 
Ss” oF 


图 1 高 能 
Figure 1 Changes in high-energy 


2.2 数据 处 理 


~~) 2 D Q] 
V A V gs 
AN a” 8 yw NS 
FS E POP E g 
v Y 


® ,9 N 
V V © L wv 
N € v g 
N © Nm 
© ANT ww 
Ff £ Ò fF LF FS 
Date 


电子 通 量 在 2006 年 12 月 19 日 -2012 年 5 月 12 日 之 间 的 变化 情况 
electron flux between December 19，2006 and May 21，2012 


由 图 1 可 知 由 于 时 间 跨 度 长 ， 数 据 量 庞大 


数据 十 分 密集 ， 高 低 起 伏 不 平 ， 直 接 用 于 分 析 研 


究 ， 很 难 直 观 的 给 出 周期 性 特征 ， 故 做 日 平均 处 理 ， 如 图 2 所 示 选 择 的 数据 是 2006 年 12 月 19 日 


到 2012 年 5 月 21 


的 通 量 是 存在 
在 2006 年 12 月 19 日 至 2008 年 12 月 1 


之 间 的 日 平均 数据 。 由 于 数据 庞大 我 们 只 选择 研究 以 “天 ” 
量 周期 且 日 平均 数据 将 平 请 掉 短 周期 时 变 特性 ， 故 采用 日 平均 数据 是 合理 的 。 该 图 显示 高 能 电子 
定 周期 的 且 周 期 不 唯一 ， 分 主 周期 和 次 周期 ， 根 据 图 像 的 分 布 特点 ， 可 以 看 1H 


为 量 级 的 电子 通 


1] 
Bj 


9 日 和 2009 年 12 H19 HÆ2012 %5 月 21 日 期 间 曲 线 波动 
较 大 ， 而 在 2008 年 12 月 19 HÆ 2009 年 12 H19 日 这 一 年 时 间 的 通 量变 化 十 分 微小 。 为 了 更 好 的 
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22.0 MeV Electron Flux 
(om?-s-sr-keV)"1 


图 2 高 能 电子 在 2006 年 12 月 19 日 -2012 年 5 月 21 日 之 间 日 平均 通 量 图 


Figure 2 The daily average flux of High Energy Electronics between December 19, 2006 and May 21, 2012 


由 于 FY-2D 卫星 的 粒子 探测 器 所 探测 到 的 高 能 电子 通 量 数据 存在 部 分 缺失 ， 为 了 便于 处 
理 ， 我 们 对 原始 数据 进行 了 插值 处 理 ， 然 后 采用 归 一 化 方法 ， 以 减少 插值 数据 对 原始 数据 的 影响 以 
及 对 整个 研究 结果 的 干扰 %。 由 图 3 可 知 数据 分 布 特点 与 图 2 所 显示 的 几乎 一 致 ， 表 明 使 用 插值 后 
的 数据 研究 周期 特性 较为 合理 。 


10 T T T T T 
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> 
让 8 上 4 
c 
Ba 
gS of 
Tg 
® 5 4 上 4 
2 4 
oN 
N 2 | 
八 e 
x 
| 0 ] 
fang 

1 


-2 1 1 1 1 
2006/12/19 2007/12/19 2008/12/19 2009/12/19 2010/12/19 2011/12/19 
Date 


图 3 归 一 化 日 平均 图 


Figure 3 Normalized daily average diagram 


3 高 能 电子 通 量 的 周期 特性 分 析 


为 估计 高 能 电子 通 量 的 周期 特性 吗 ， 本 文选 取 2006 年 12 月 19 日 -2012 年 5 月 21 日 的 日 平 
均 数 据 进 行 分 析 。 分 别 采 用 最 大 信 谱 法 和 Levinson-Durbin、Burg 算 法 进行 计算 ， 对 比分 析 
两 种 方法 的 准确 性 。 
3. 1 AR 模型 阶 数 选择 与 确定 

由 2.2.1 可 知 ， 在 AR 模型 参数 的 选择 时 ， 我 们 一 般 采 用 FPE 和 AIC 两 种 准则 来 判断 。 
选择 正确 的 AR 阶 数 是 最 大 科 谱 估 计 的 关键 之 一 ， 其 值 的 大 小 能 够 影响 我 们 判断 箭 谱 图 的 波 
峰 ， 继 而 影响 周期 的 研究 结果 。 

根据 FPE, AIC 准则 的 计算 公式 ， 基 于 matlab 来 实现 最 佳 阶 数 的 选择 [11] ， 如 图 4 所 示 ， 
4( 左 )、 图 4( 右 ) 分 别 展示 的 是 FPE 准则 和 AIC 准 则 关于 阶 数 和 参数 估计 的 关系 ; 


202011.00015v1 


chinaXiv 


FPE parameter estimation 


一 * 一 FPE 
~~ Minimum of FPE 


ChinaXiv 合 作 期 刊 


x104 


—+— AIC 
3.3027 Minimum of AIC 


7 8 9 10 11 12 


BAE 4) PA 


Order of AR model 


图 4 FPE 准则 (AE) 和 AIC 准则 CH) 


Figure 4 FPE guidelines (left) and AIC guidelines (right) 


PHBE. BA 


标 分 别 为 参数 估 值 和 所 对 应 的 阶 数 。 用 时 间 序 列 的 最 小 


二 乘法 和 FPR、AIC 准 则 对 AR 模型 进行 参数 估计 ， 再 作 折线 图 得 到 对 应 的 FPE 和 AIC 准 则 函数 关系 。 
由 图 4 ( 左 ) 、 图 4 CA) 两 图 对 比 可 以 看 出 AR 模型 阶 数 为 10 时 ，FPE 和 AIC 函 数 所 对 应 的 参数 佑 值 
均 为 最 小 值 ， 分 别 为 1277145. 757 和 33018. 72179, 根据 上 面 所 说 的 确定 标准 ， 运 算 结 果 符 合 
(11) 。 由 此 确定 该 AR 模型 的 最 佳 阶 数 为 10。 所 确定 的 最 佳 阶 数 均 适 用 于 下 面 
AY fee Ae ih tt ALevinson-Durbin, Burgtiyž. 


EX2. 3 公式 


3.2 高 能 电子 通 量 周期 特性 分 析 


3.2.1 RAMI it 
根据 4.1 WARES EEA, EAT ay Tee). A AP i (He FE FH ec tc 
KRENEK, WE 5 所 示 ， 阶 数 为 10， 高 能 电子 通 量 周 期 分 主 、 次 周期 na， 根据 功率 谱 
概率 分 布 可 知 ， 分 布 图 存在 有 多 个 峰值 ， 本 文选 择 研究 的 是 最 高 波峰 和 次 高 波峰 。 最 高 波峰 频率 为 
0.0625d-! 时 ， 对 应 为 主 周 期 ， 功 率 谱 最 大 ， 值 为 14.16， 对 应 周期 为 13.87d， 次 波峰 频率 
0.1641d' 时， 对 应 为 次 周期 ， 功 率 谱 处 于 第 二 大 值 ， 值 为 13.71， 对 应 周期 为 27.8d。 此 研 
究 结 果 与 相关 文献 叫 利 用 小 波 分 析 得 出 的 高 能 电子 通 量 的 27d 周 期 和 13d 周 期 结果 几乎 一 
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外 子 通 量 周期 具有 一 定 准确 性 。 
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图 5 eK MTR 


Figure 5 power spectrum diagram of Maximum entropy spectrum method 


算 结果 大 概 存 如 


E0. 03-0. 08d 的 误差 。 


此 外 图 中 可 


也 能 判断 出 其 不 稳定 性 9。 


EE 可行， 也 考虑 了 诸多 因素 ， 但 是 由 于 对 数据 进行 了 插值 处 理 ， 经 初步 分 析 ， 计 


以 明显 看 出 功率 谱 不 止 一 个 波峰 ， 表 明 周期 不 唯一 ， 且 根据 图 像 波峰 变换 的 趋势 
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3.2.2 Burg 算法 和 levinson 算 法 

上 文 已 经 确定 了 最 佳 阶 数 为 10， 采 用 日 平均 处 理 后 的 数据 进行 周期 特性 分 析 。 首 先 按 照 最 
佳 阶 数 10 绘 出 Burg, Levinson 算法 的 功率 谱 图 像 ， 与 图 5 对 比 ， 图 6 的 左 、 右 两 图 像 基 本 一 致 ， 再 次 
证 明 阶 数 10 为 最 佳 阶 数 。 而 观察 左 、 右 图 可 知 主 周期 频率 为 0.011d1 ， 次 周期 频率 是 0.128d- | 所 算 
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图 6 ”Burg 算 法 (Æ) ~ Levinson 算法 〈 右 ) 功率 谱 图 


Figure 6 power spectrum diagram of Burg algorithm (left) and Levinson algorithm (right) 


根据 上 两 图 对 比 可 知 ， 高 能 电子 通 量 确实 存在 周期 性 ， 且 有 多 个 周期 ， 周 期 分 为 主 、 次 周 
肯 ， 分 别 为 13. 87d 和 27. 8d。 与 此 同时 ， 对 比 了 最 大 箭 谱 法 和 Burg. Levinson 算法 的 分 析 
寺 果 ， 结 果 显 示 最 大 峭 谱 法 的 分 析 结 果 更 接近 真实 值 ， 更 符合 实际 情况 。 

相关 文献 表明 ts 可 ， 电 子 通 量 的 13 天 和 27 天 周期 与 太阳 活动 密切 相关 。 比 如 在 太阳 活 
低 年 ， 横 越 赤道 的 冕 洞 发 展会 产生 高 速 太 阳 风 ， 导 致 外 辐射 中 电子 通 量 发 生变 化 ， 随 着 太阳 的 自 
转 ， 延 伸 型 的 冕 洞 每 27 天 旋转 一 周 ， 从 而 电子 通 量变 化 也 具有 了 27 天 的 周期 性 。 而 电子 通 量 的 
13 天 周期 也 是 伴随 着 太阳 风速 的 变化 而 产生 的 ， 在 27 天 周期 中 主要 有 3 个 相 期 : 快速 上 升 、 峰 
值 期 和 下 降 期 ， 这 种 对 称 变化 过 程 ， 即 形成 了 13 天 周期 ， 所 以 一 般 在 电子 通 量 的 27 天 周期 中 
存在 13 天 周期 的 变化 成 分 。 


4 结论 


NI 
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本 文 基于 风云 二 号 D EARME TT EHS AOE, BB OTE ZO a PEMA TST 
究 ， 得 出 以 下 几 个 结论 : 
(1) 最 大 炉 谱 估计 表明 高 能 电子 通 量 的 主 周期 为 13. 87d 和 27. 8d， 而 且 高 能 电子 通 量 的 周 
期 不 唯一 ， 不 稳定 。 

(2) JOT Levinson-Durbin、Burg 算法 关于 高 能 电子 周期 特性 研究 的 比较 ， 
可 以 判断 出 最 大 粒 谱 估计 更 适合 于 高 能 电子 周期 的 研究 ”“， 结 果 更 接近 于 真实 情况 。 

(3) 太阳 活动 会 影响 高 能 电子 的 周期 ， 随 着 太阳 活动 的 变化 ， 周 期 也 将 发 生 相 应 变化 。 通 
过 分 析 高 能 电子 通 量 的 周期 变化 可 以 侧面 预测 太阳 活动 的 变化 趋势 ， 为 日 后 研究 太阳 活动 提供 有 
价值 的 参考 。 
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Jiading 201800, China;) 
Abstract: Based on maximum entropy spectral estimation method, this paper analyzes the 

periodic characteristics of the high-energy electron flux 22.0MeV which has been observed by 
FY-2D meteorological satellite. The algorithm is based on the autocorrelation model (AR), a period 
of 13.87d and 27.8d is deduced by analyzing the power spectrum. According to the FPE and AIC 
criteria, the parameters of the AR model can be determined. Levinson-Durbin algorithm and Burg 
algorithm are compared with the maximum entropy spectral estimation method, it is found that 
maximum entropy spectral estimation has advantages in the study of periodic characteristics. The 
results are very important for the study of the spatial distribution of high-energy electrons in 
geosynchronous orbit,to forecast high energy electron flux enhancement events,or early warning 


of deep charging events. 


Key words:High-energy electrons; Maximum entropy spectral estimation;The periodic 


characteristics; power spectrum estimation; 
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